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Abstract 

We study the effect of the electron-electron interaction on the transport of 
spin polarized currents in metals and doped semiconductors in the diffusive 
regime. In addition to well-known screening effects, we identify two additional 
effects, which depend on many-body correlations and exchange and reduce the 
spin diffusion constant. The first is the "spin Coulomb drag" - an intrinsic 
friction mechanism which operates whenever the average velocities of up- 
spin and down-spin electrons differ. The second arises from the decrease 
in the longitudinal spin stiffness of an interacting electron gas relative to a 
noninteracting one. Both effects are studied in detail for both degenerate and 
non-degenerate carriers in metals and semiconductors, and various limiting 
cases are worked out analytically. The behavior of the spin diffusion constant 
at and below a ferromagnetic transition temperature is also discussed. 

72.25.Dc, 72.25.-b, 75.40.Gb 
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I. INTRODUCTION 



The theory of spin-polarized transport is attracting the attention of the physics com- 
munity both for its potential applications to the emerging field of "spintronics" [|J and in 
relation to the intriguing field of quantum computation 0. In this context Kikkawa et al. 
( OS) generated much excitement by demonstrating optical injection and subsequent con- 
trol of packets of spin polarization in the conduction band of n-doped GaAs. Such packets 
could in principle be used to transport information between separate regions of a system. 
Unlike ordinary electron-hole packets, whose mobility is limited by strong scattering in the 
hole component, these unipolar electronic spin packets are both long-lived (with lifetime 
t s ~ 10 ns) and highly mobile M. 

From the theoretical point of view, most recent work on spin transport phenomena is 
based on an independent electron model, where, in general, the coupling between different 
spin channels is completely neglected. Flatte and Byers |5|] have recently discussed the 
behavior of spin packets in semiconductors in the framework of the Hartree approximation, 
where the Coulomb interaction is taken into account only through the imposition of a charge 
neutrality constraint. This constraint establishes an electrostatic coupling between different 
spin channels. They neglect, however, all further many-body effects. 

In this paper we go beyond the treatment by Flatte and Byers by considering more subtle 
correlation (spin-drag) and exchange (spin-softening) effects. Our objective is the derivation 
of drift-diffusion equations for spin packets in a full many-body context, i.e., allowing for 
correlation between different spin channels. These equations contain the spin-packet mobility 
and the diffusion constant as key parameters. We show that both exchange interactions and 
the drag effect between different spin populations ||, concur in reducing the value of the 
diffusion constant by a sizable amount, leaving at the same time the mobility basically 
unaffected. The microscopic quantities responsible for this effect are the longitudinal spin 
stiffness (the second derivative of the free energy with respect to magnetization, not to be 
confused with the transverse spin stiffness of the ferromagnetic state) and the spin-drag 
transresistivity. We discuss in detail the behavior of these quantities in various regimes and 
show how important the Coulomb effects are when dealing with doped semiconductors (see 
Sec. HD. 



This paper is organized as follows: 

In Section II we review the basic ideas underlying the spin-resolved drift diffusion equa- 
tions and make use of the Landau transport equation to elucidate the structure of the 
homogeneous spin resistivity matrix. In particular, we show that the off-diagonal element 
of the spin-resistivity matrix (the spin trans-resistivity) is almost exclusively controlled by 
the Coulomb interaction: the contribution from spin-flip scattering, while finite, is utterly 
negligible for short-range scatterers. 

In Section III we present the calculation of the Coulomb contribution to the spin trans- 
resistivity. The theory of Ref. || is extended in various directions. First we study the 
spin drag effect as a function of temperature going from the degenerate regime (which is 
appropriate for ordinary metals) to the nondegenerate regime, which is appropriate for low- 
density/high -temperature doped semiconductors. Then, we calculate the spin drag in the 
"mixed" case in which one spin component is degenerate while the other is nondegenerate: 
this is relevant to situations in which a strong spin polarization exists. 
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Section IV is devoted to a description of the behavior of the longitudinal spin stiffness 
of the homogeneous electron g function of density and temperature. 

In Section V we present a detailed derivation of the drift-diffusion equation for a macro- 
scopic spin packet, and give explicit expressions for the mobility and diffusion constant in 
terms of microscopic quantities such as the spin stiffness and the spin transresistivity. 

Section VI suggests experiments aimed at directly measuring the spin transresistivity 
and the spin diffusion constant in metals or semiconductors. 

Electron gas theory predicts that, at sufficiently low temperatures and densities, the 
homogeneous electron gas undergoes a ferromagnetic transition |7]||. Such a transition could 
be in principle observed in a doped semiconductor . Section VII examines the behavior of 
spin diffusion constant in the paramagnetic phase and how the occurrence of ferromagnetism 
as a second order phase transition would affect it. Our calculations show that the diffusion 
constant vanishes at the transition temperature and increases as the system becomes fully 
spin polarized with decreasing temperature. 



II. GENERAL THEORY I 

A. The drift-diffusion equation 

The theory of diffusive transport in metals and semiconductors is based on the assump- 
tion of a local linear relationship between the current densities of up and down spin electrons 
and the gradient of the local electro-chemical potentials: 

- eJ a (f) = - £ a a/3 (f)V^(f). (1) 

P 

Here e is the absolute value of the electron charge, J a {r) (a =f or 1) are number currents 
( H) and the electro- chemical potentials ip a ( r ) are given by the electrostatic potential <f)(r) 
plus the local chemical potential, which can be spin-dependent, 

=*<?)- ~—Q^ ■ (2) 

Here f(n^,n^,T) is the free energy per unit volume of a homogeneous interacting electron 



gas evaluated at the local spin densities n a (f) and uniform temperature T JT0[; a a/3 = 
lim^o cr a/3 (g = 0,u) is the homogeneous conductivity matrix of the electron gas, whose 
structure, in the presence of interactions, will be elucidated below. 
Substituting (JJ) in (P and writing 

^ df(n h n b T) d 2 /(n T , ra i; T) ^ 



we obtain 



ef a {r) = E (Mr)V0(r) - eD^(r)Vn p (f)) (4) 



P 

where the diffusion matrix D a p is given by 
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e 2 D a/3 = ^2a a ^S 1/3 (5) 

7 

and 

_ 9 2 /(n T ,n i) T) 
a/? dn a drip 

is the static longitudinal spin-stiffness matrix - the inverse of the spin susceptibility matrix 
Xa/3- Eq. (^) is the well known Einstein relation between the mobility and the diffusion 
constant, generalized here to the case of spin polarized transport. The first term of Eq. (|J) 
is the drift current associate with the electrostatic field, the second is the diffusion current 
associated with the gradient of the electronic densities. These two terms cancel out exactly in 
a situation of local equilibrium, due to the Einstein relation and the fact that = S a p. 

On a formal level the main effect of the Coulomb interaction is the appearance of non- 
vanishing off- diagonal elements of the conductivity and spin- stiffness matrices, ^ 
implies that an electric field acting only on the up-spin electrons must necessarily drag 
along a current of down-spin electrons. Conversely, a current of up-spin electrons Jf ^ 
flowing against a background of stationary down-spin electrons (Jj_ = 0) will necessarily 
induce a gradient of spin-down electro-chemical potential = —Vipi = ep^J^, where p Q/3 
is the resistivity matrix, inverse to o~ a p. We shall see later how one can make use of these 
effects for a direct experimental measure of the spin Coulomb drag. 

The other noteworthy feature is Sfj, 7^ 0: this means that the chemical potential of up- 
spins df{n^ n^T)/ dn^ is a function of both up and down spin densities. Thus, a disturbance 
acting on one of the two spin populations will affect the other through Coulomb correlation. 



B. Structure of the resistivity matrix 

Although the homogeneous resistivity matrix can be calculated from first-principle Kubo 
formulas and/or transport equations its general structure (including the off-diagonal terms, 
due mostly to the Coulomb interaction) is best understood at the phenomenological level. 
We first present the phenomenological picture and then justify it from a more formal con- 
sideration of the Landau transport equation. 

Let E^(t) and Ei(t) be uniform effective electric fields, E a = — V^ Q , that couple to up- 
and down-spins respectively. We restrict ourselves to the linear response regime. If v a is the 
velocity of the center of mass of electrons of spin a, and N a the number of such electrons, 
then the equation of motion for v a has the form 

-> -* in* tti* 
m*N a v a = -eN a E a + NJ a + — N & v & , (7) 

where m* is the effective mass of the carriers and a = —a. 
Let us examine the meaning of the various terms in ([7|). 

The first term on the right hand side is the net force exerted by the electric field on 
spin-a electrons. 

The second term, F^ & , is the Coulomb force exerted by spins of the opposite orientation 
a on spin-a electrons. Notice that the net force exerted by spins of the same orientation 
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vanishes by virtue of Newton's third law. For exactly the same reason we must have = 
—F^ a , and by Galilean invariance this force can only depend on the relative velocity of the 
two components. Hence, in the linear approximation, we can write 

Fg i = - 1 mTN a 1 ^(v a -v ei ), (8) 

where n = + nj_ is the total density. Eq. (§) defines the spin drag coefficient 7. 

The third term on the right hand side of (0) combines two distinct physical effects. One 
is the net force exerted on spin-a electrons by electron-impurity collisions that do not flip 
the spin of the incoming electrons. The other is the rate at which momentum is lost to 
the a component as a result of electron-impurity collisions that flip the spin from a to a. 
Accordingly, the momentum relaxation rate r" 1 is written as the sum of non-spin-flip and 
spin-flip contributions: r" 1 = T~t a + T~} a . Spin-flip times have been found to be very long 
both in metals (r s / )Q ~ 10 -1 ns and in semiconductors (r s / )Q ~ 10ns 0,|J), whereas 
non-spin-flip times are usually much shorter (of the order — 10~ 4 ns). Thus, for most 
practical purposes, r^ 1 ~ T nf,a- 

Finally, the last term on the right hand side of (0) represents the rate at which momentum 
is fed into the a component by electron-impurity collisions that flip the spin from a to a. 
Because electrons emerging from collisions have a randomized momentum distribution with 
nearly zero average, we expect (t^) _1 to be smaller than even the already small spin-flip rate 
T~f a . Indeed, we shall see in the next section that (t^,) -1 vanishes for short-range scatterers 
in the Born approximation. 

Fourier-transforming the equation of motion ([7p for the current density j a (oj) = 
—en n v n (ui) we find 



iujj a (uj) 




(9) 



The resistivity matrix p aa i is defined as the coefficient of proportionality between the electric 
field and the current, i.e. E a = J2 a ' Paa'ja'- A quick comparison between this definition and 
Eq. (||) shows that the complete form of the resistivity matrix p a p is 

/ -lu^f- + + ~5~t - ^7 \ 

P=\ »^_»il 7 _ iu) m!L + + TUm^y I ' ^ 1U > 

y n^e A T^ ne z ' e n l n l e T l n l ne / 

Notice that this matrix is symmetric, due to the relation \/n]T t \ = l/n^ri, which will be 
proved in the next section. 

Due to the extreme smallness of the spin-flip rates l/r' a , the off-diagonal resistivity (also 
referred to as "spin-transresistivity" ) is controlled almost entirely by the Coulomb interaction 
term, i.e. we can safely assume 

771* 

PU = 2^ ( H ) 

ne 2 

and 7 is directly proportional to the spin trans-resistivity. But then, Galilean invariance and 
Newton's third law demand that the same appear also as a correction to the ordinary 
diagonal resistivity. This effect is quite distinct from the "trivial" renormalizations due to 
the electronic screening of these interactions. 
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C. Derivation from Landau transport equation 



In this section we provide a microscopic justification of the phenomenological Eq. (^) 
and give explicit expressions for the electron-impurity relaxation rates. 

We start from the linearized transport equation for the quasiparticle distribution function 
in Landau theory of Fermi liquid |12|]. This is simply the classical Boltzmann equation for 
quasiparticles in a self-consistent field described by Landau parameters. In the homogeneous 
case the self-consistent field coincides with the classical electric field and the distribution 
function f a (k,t) obeys the well-known kinetic equation 

9fa(k,t) c ^ dfWjk) f df a (k,t) ^ 



dt dk \ dt 



coll 



where f£\k) is the equilibrium distribution function. 

The collision term (df a (k,t)/dtj ^ has contributions from the Coulomb interaction as 
well as spin-flip and non spin-flip electron-impurity interactions. The various contributions 



are listed below 



d f a (k,t ) 

dt 



Coulomb collisions 

c 

= - J2 wC {ka,p(3]k'a\0(3')5j: +i ;j: l+f 5 a+ p^ + pi 

coll pk'f,l3a'l3' 

{Uk)U{m - fAk'W - fptf)] - [1 - /<*(£)][! - Uip)]fAk')U'i0)} 

8(eka + — Ck',/3' — e p'P')i (13) 

where W c (ka,pj3; k'a' ,p* (3') is the probability of the Coulomb scattering process 
ka,pP — > k'a', ftp' and e^ a is the energy of a particle of momentum k and spin a 
relative to the chemical potential. The conservations of momentum, energy , and spin 
are explicitly displayed. 

Non-spin-flip electron-impurity collisions 

nf 

= -^2W nf (ka,k'a)[f a (k) - f a (k')]S(e ka - £k , a ), (14) 



df a (k,t) 



ot 



coll 



where W n * (ka, k'a) is the probability of the non-spin flip (nf) scattering process ka 
k'a. 



Spin-flip electron-impurity collisions 

= -J2 Wsf tfa,k'a)[fa(k) - fa(k')]8(e ka - e# a ), (15) 



df a (k,t) 



coll 



ot 

where W s f (ka } k'a) is the probability of the spin flip (sf) scattering process ka — > k'a 
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In order to obtain a closed equation of motion for the currents, such as Eq. @, we must 
multiply both sides of Eq. (|E3) by —ek/m* , sum over k, and then express the integrated 
collision term 




df a (k,t ) 
dt 



(16) 

coll 

back in terms of the currents. Of course, this cannot be done rigorously, but for an isotropic 
system slightly perturbed from equilibrium one can assume Jl3| that the distribution function 
of the state with currents j a = —en a v a is given by 

failt) = fi°\e ka ) - df l {€ka) v a (t) ■ k. (17) 

Substituting this into Eqs. (|13|), (HI), an d (HH), an d linearizing with respect to the currents 
wherever needed, we arrive, after tedious but straightforward manipulations, at the desired 
equation of motion (||), with the following expressions for the various relaxation times: 

7 = m *jV a n„ JS 2dk B T W ^ ^' k ' a ' V ^) 6 k+p,k'+p' 6 ( e ^ + ~ £k>,a - e P >a) 

kpk'p' 

fL \e ka )fL \e pa )fi \-e k , a )fi°\-e p , a ), (18) 
where d is the number of spatial dimensions; 

— = - E df ^ {tka) — E W n f(ka, k'a)(l - k ■ k')5(e kQ - e Wa ), (19) 

^ ^ — * — * 

where k and k! are unit vectors in the directions of k and k!\ 

I ^ dfjp (e ka ) e ^ waf ^ %a)5(e ka - eya ), (20) 



7~sf,a j* d^kct Tl a d 

\ = ~ E df l (eka) E W sf (ka, k'a)k ■ k'5(e ka - e Vs ). (21) 

Notice that l/r s f >a and l/r„ arise, respectively, from the first and the second term on 
the right side of Eq. (|T5|). The key difference between these two relaxation rates is that 
the expression for the latter involves an angular average of the scattering probability with 
weight factor k ■ k = cos(8). This average vanishes in the Born approximation for short- 
range scatterers, since the scattering probability becomes isotropic (independent of k and 
k') in this special case. Quite generally, one can expect l/r' a to be much smaller than l/r s / )Q 
in agreement with the qualitative arguments given in the previous section. This means 
that the spin transresistivity is almost entirely a Coulomb interaction effect, and therefore 
its measurement can shed light on the nature of the Coulomb correlation between up- and 
down-spin electrons. This is one of the main points we wanted to make in this section. 

Finally, notice that Eq. (|21|) implies the identity 

A = (22) 

n i T i n T r T 

which guarantees the symmetry of the resistivity matrix Eq. (|l0"D . 
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III. CALCULATION OF THE SPIN TRANSRESISTIVITY 



The theory of the spin transresistivity has been worked out in ||. This theory closely 
parallels the theory of the ordinary Coulomb drag between parallel two-dimensional electron 
or hole-gas layers but differs in some important details, as the fact that electrons of 
opposite spin interact with the same set of impurities, so that certain electron-impurity 
terms which appear in the Kubo formulation of the transresistivity do not vanish upon 
disorder averaging. Fortunately, it turned out that these terms cancel out exactly at low 
frequency (uj « E F ) and to leading order in the electron-electron and electron-impurity 
interactions ||. 

In this section we first review for completeness the derivation of Ref . , and then present 
an RPA calculation of pn at finite temperature At variance with the calculation of 



Ref. || we present our results not only in the low temperature limit (which is relevant to 
metals and where 7 ~ T 2 ), but also in the non-degenerate k B T » E F and quasi-degenerate 
k B T ~ Ep regimes (which are relevant to doped semiconductors), where k B is the Boltzmann 
constant and Ep = h 2 (37i 2 n) 2 ^ 3 / 2m* is the Fermi energy ||16|| . 

We start from the Kubo formula |T7| for the uniform conductivity matrix 



(?a,a>{u) = r n a 6 a>a/ H , (23) 

iuj m* \ m* ' 



where ((A; B))^ represents, as usual |]T7[, the retarded response function for the expectation 
value of A under the action of a field that couples linearly to B. We assume to be in the 
diffusive, weak scattering regime characterized by h/r D « k B T, where t d = {n^/njr^ + 
(nJnjT^. Because of the "high" temperature, weak localization effects are negligible. In this 
regime the resistivity is essentially independent of frequency for uj « k B T /Ti. It is therefore 
legitimate to take the limit of weak electron-impurity and electron-electron scattering before 
taking the limit of uj — > |Tj|. When the limits are carried out in this order the P a 's are 
almost constants of motion and therefore the second term in the square bracket of Eq. ( |23| ) 
is a small correction to the first. Inverting Eq. (|23D to first order in ((P a ; P a >))u an d selecting 
the t| matrix element we obtain 

= ^«M>k. (24) 

e z n^ni 

It is convenient to recast this equation in a form that emphasizes the importance of the non 
conservation of Pj and P|. To this end we make use twice of the general equation of motion 

((A ] B)) u = -{{[A 1 B\)+i((A;B)) u ) 1 (25) 

where A = —i[A, H] is the time derivative of the operator A, and (..) denotes the thermal 
average. Thus, Eq. fl2"4j) can be rewritten as 

(w) = i (<M>k±i<!Ml>. (26) 
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The commutator term controls the high frequency behavior of p-\i{ui) and can be expressed 
in terms of ground-state properties ||19|| . This term however gives a purely imaginary contri- 
bution to the trans-resistivity Our present interest is in the real part of the trans-resistivity, 
which is controlled by the imaginary part of the force-force response function. 
The force operator is given by 

^ ^ 

P ^ = -yJ2 WqP&P-fa ~ y £ Q V 7~ V^-tf* > ( 27 ) 

where v q = 47re 2 /q 2 e is the Fourier transform of the Coulomb interaction, e the dielectric 
constant of the material, v e ~ % is the Fourier transform of the electron-impurity interaction, 
Pq a is the electronic spin density fluctuation operator, is the Fourier transform of the 
impurity density (a number), and V is the volume of the system. 

As explicitly shown in ||, the contribution of correlated impurity scattering to the trans- 
resistivity, in the low frequency Tiu « Ep limit and to leading order in the electron-electron 
and electron-impurity interactions, vanishes since the Coulomb-impurity term exactly cancels 
the impurity-impurity contribution. Thus, in this limit, the real part of the spin trans- 
resistivity takes the form 



1 ^ 



Re PnH = ——a, ,t/2 £ 



n^n^e 2 ujV 2 3 



-VqVq' 



i™{{p-<nPq~i;p<^p-<n))u- (28) 

We proved the cancellation of correlated impurity scattering effects within the frame 
of the Drude-Boltzmann theory defined by Eq. ( p4j) which is the result of interchanging 
the natural order of the u — > limit and the weak scattering limit. This approach is 
only adequate in the "classical" regime ksT » Ti/td- A more sophisticated treatment 
of quantum effects in correlated impurity scattering |2(J suggests that the spin drag would 
be even larger than predicted by the present theory at temperatures ksT << H/td- The 
temperature range in which these quantum corrections are important shrinks to zero in the 
limit of weak impurity scattering. 

We have calculated the four point response function 
X4 P (<f, <f, to) = {{p-q\Pqi, Pq> tP-cf i))ui at finite temperature in a generalized Random Phase 
Approximation (RPA). Because of its infinite range, the Coulomb interaction must be treated 
to infinite order, even when weak. The sum of the RPA diagrams M has been evaluated by 



standard methods [21| with the following result: 



1 1 n 2 (p-P^ 



n^n^e 2 V ~? 3 9 uj 
-oo 7i \e~W - l)(e-^— ') - 1) 



(29) 



Here (3 = 1/ksT, Xaa'iQ^) * s the imaginary part of the RPA spin-resolved density- density 
response function, which is related to the noninteracting response function xoa{c[,u) as 
follows 

[X~ l (q, U))] aa > = [XOaY 1 ^ ^)<W ~ V q . (30) 
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It is possible to show by simple but tedious algebraic calculations that Eq. (|29|) for the 
spin trans- resistivity p^(u>,T) reduces, in the case of finite temperature and uj = 0, to 
the well known result of memory function and diagrammatic theories for the conventional 



Coulomb drag [g|], |3 



HePU ^ 1] -^^v\ J V «2Jo V| e (g,^)| 2 sinh 2 (^'/2)' { ' 

Furthermore, for T = and u ^ 0, the RPA is equivalent to the decoupling approxima- 
tion for the four-point response function used in |23j] to calculate the dynamical exchange- 
correlation kernel. In this limit, the real part of the spin transresistivity takes the form 

Rep n (u, 0) = \ v l 

o d ^-\ X !^{q,u')-x!{ i {-q,u - uj') - X "i(q, tj')^{-q, uj - uj')] (32) 

Thus our calculation demonstrates that those two approximations, quite different at a first 
sight, are simply RPAs performed in different limits. 

From now on, for simplicity of notation, we will refer to Rep^i simply as p^. 



A. Numerical evaluation 



To calculate at finite temperature, we have used in Eq. fl2"9|) the temperature de- 
pendent expression for the three-dimensional noninteracting spin-resolved density-density 
response function 



xLi^^T) 



1 



1 



167T qa* 3 Ry 



uj 



1 . l + e 
In 



PRy 



l + e 



j_/ Sa>-e g V 
^\ 2 J 



(33) 



where a* is the effective Bohr radius, q = qa*, uj = hu/Ry, Ry = e 2 /2a* is the effective 
Rydberg, £ Q is the chemical potential for the a spin population and e q = h 2 q 2 /2m*. Eq. 
follows directly from the definition 



k 



(34) 



where n^ a = l/[exp(/3(e g — £ a )) + 1] is the average number of a-spin electrons with energy 



e q . 



Fig.g shows \pii\ as a function of temperature and density pSJ . The data are calculated in 
the paramagnetic phase and for semiconductor parameters (GaAs), i.e. m* = 0.067, e = 12 
and carrier density ri\ — 1.5 x 10 16 cm -3 , ri2 = 1.5 x 10 17 cm" 3 , and n.3 = 1.5 x 10 18 cm -3 , p^ 
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peaks at about the Fermi temperature Tp, underlying the crossing between the degenerate 
and the non degenerate regimes. As can be seen, p-j-j is strongly enhanced as the density 
decreases, mainly due to the prefactor dependence ~ 1/n 2 . In fact its maximum increases of 
almost two orders of magnitude, passing from 0.3 millifi x cm for n = n\ to 14 millifi x cm 
for n = 713. In the calculations of the following sections, we will mainly focus on the density 
value n — n,2, corresponding to a Fermi temperature Tp = 178K. The inset of Fig. [I] 
presents for this density value the comparison between pjj_ and its non-degenerate analytical 
approximation (dashed lines) discussed in details in Sec. |III C| . 

We now turn to a quantitative assessment of the relevancy of the spin Coulomb drag. 
First of all it is necessary to underline that the spin drag is an intrinsic effect of spin- 
polarized transport: that is, while impurity scattering could in principle be suppressed in a 
perfect crystal, the spin Coulomb drag will always be present, even in the purest sample, and 
dominate over phonon scattering at sufficiently low temperature. However, since available 
samples are usually far from perfection, it is reasonable to ask how the spin transresistivity 
compares to the more familiar Drude resistivity. In metals, as we shall show in detail in 
the next section, one finds, at most, pn ~ 10 _2 /if2 x cm so that p^J Pd is of the order 
of few percents. The situation is very different for semiconductors: since both the Fermi 
temperature (at which pjj_ peaks) and the carrier density are considerably lower than in 
metals, pjj_ can become comparable and even greater than pp. This strong variation depends 
on specific semiconductor characteristics, such as the effective mass, mobility, and density 
of the carriers. 

In Fig. H we show the effect of the carrier mobility p on the ratio Ph/pd- We plot Pu/Pd 
in respect to temperature for an n-doped semiconductor as GaAs (m* = 0.067m e , e = 12 
and 11 = 1.5 x 10 17 cm -3 ; upper panel) and for a p-doped semiconductor, as (Ga,Mn)As 
(m* = 0.5m e , e = 12, n = 1.2 x 10 19 cm -3 ; lower panel). Each curve corresponds to a 
different mobility value as reported in the figure caption. The values increase from "A" to 
"D". In particular the value p = 3 x 10 3 cm 2 /Vs (labels as "C"), corresponds to the value 
measured for a spin packet in 0]. As can be seen, changing the material it is possible to 
increase the ratio p^J Pd by an order of magnitude, to the point that the spin transresistivity 
can become greater than the Drude resistivity. 

In Fig. |]we show the dependence on the ratio p^J pp> on the carrier density for GaAs. The 
results are presented for two different temperatures (T = dashed lines and T = 300i^, 
solid line) and two different values of the mobility (as labelled in the figure). Each curve 
peaks at a density such that T F = T. 

In conclusion, our calculations demonstrate that in semiconductors, at T « Tp, the spin 
Coulomb drag must definitely be considered an important contribution to the resistivity for 
spin polarized currents. 

B. The degenerate limit 

Let us now focus on the low-temperature (&bT << Ep) and low-frequency (tku « Ep) 
regime. This is the only regime of practical importance in ordinary high- density metals. 
The spin-Coulomb drag coefficient is controlled by a subset of all the processes that lead to 
the finite lifetime of a quasi-particle at the Fermi surface, namely, the processes in which the 
quasi-particle in question exchanges momentum with an electron of opposite spin (scattering 
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processes between parallel spin electrons do not cause relaxation of the spin-current). Since 
the inverse quasi-particle life-time at the Fermi surface is known in Fermi liquid theory to 
scale as (k B T / Ep) 2 we expect the same scaling to hold for the spin drag coefficient at low 
temperature. 

This prediction is confirmed by the detailed calculation as follows. At low temperature, 
the exponential factors in Eq. (p9|) restrict the region of integration to uj ~ k B T/h. The 
low-frequency form of the density fluctuation spectra Xaa'(.Qy u ) * s a linear function of to. 
In the limit of vanishing impurity concentration XoaiQ-,^) is simply the Lindhard function, 
whose imaginary part, at low frequency, is given by Xo a (& u — ► 0) = — (m* 2 / Ak){uj / q) and 
whose real part can be approximated by its value at uj = 0. Making use of these limiting 
forms, the calculation of pfj_ can be carried in an essentially analytical fashion. The result is 



D , _ ha*4n 2 (k B T) 2 + h 2 u 2 

R ePn(^T) = -— _ 

r 2k F a* dq 



2Air 3 nin^ Jo q 2 \e(q/a* 



(35) 



where kp = min(kp^, kp^), with kp a the a spin population Fermi wave- vector, 
and e(q, uj) = 1 — f^Xot (<?> u ) ~ v qXoi(Qi UJ ) i s the RPA dielectric function. Eq. (j3~5| ) shows 
that, in the absence of impurities, ph(uj,T) is proportional to uj 2 for k B T <C Tiuj and to T 2 
for huj <C k B T. 

Modifications in the form of Xoa{q,w) due to the presence of impurities can be taken 



into account through Mermin's approximation scheme fL3|| . These modifications amount 
to replacing ujjqvp by uj/Dq 2 (D = vpr/3 being the diffusion constant) for uj < 1/r and 
q < 1/vpT where vp is the Fermi velocity and r is the electron-impurity mean scattering 
time. The uj and T dependencies of Eq. ([35]) are not affected. 

Writing explicitly in Eq. ( |5o] ) the dependence on r sa (where r sa = (47m Q a* 3 /3) _1//3 is the 
value of r s for spin a) one can also see that ph(uj,T) ~ rf+rfi, so that |p|J will strongly 
increase with decreasing electron density. In Fig. |we plot |ptj.(u; = 0, T)| as a function of 
the temperature, for n-f = ni and in the density range 1 < r s < 7. The Figure shows that, 
for metallic densities correspondent to r s ~ 5 and temperatures of the order of 40 — QOK (at 
which for example experiments on spin relaxation time using spin polarized currents have 
been performed |]TT||) , the spin trans-resistivity is appreciable [\p^{uj = 0,T)| ~ O.Olpficm). 



The corresponding Coulomb scattering time is 7 _1 from Eq. For r s = 5 we obtain 

7 _1 w 10 _13 s and S s /vp ~ 10 _10 s: 7 _1 is indeed several orders of magnitude smaller than 
the spin-flip time. This demonstrates that neglecting spin-flip processes is indeed a good 
approximation for this kind of metals. 



C. The non-degenerate limit 

The non-degenerate limit is characterized by T >> Tp. First of all we calculate 
the non-degenerate limit of the non-interacting temperature-dependent spectral function 
Eq. (^): starting from the definition Eq. (0), we use the classical expression for the fugac- 
ity exp(/3£ Q ) = n a ■ 87r 3 (/3/2m*7r) 3//2 , and obtain 
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„ / r r i \ V2KPm*n a 3e q 8h 2 u 2 . ifiw. 
XoaKQ, w; T) = — exp( — — ) exp( -— ) smh( — -). (36) 



hq v 4 7 v 4e 



In addition, in order to calculate the non-degenerate limit pjj_, we have used the classical 
limit for the dielectric constant e(q,u) = 1 + (A7ce 2 /eq 2 )(n/k B T). The final result is 



ie 2 ^m* f°° xexp{—x) 



le 2 \/m* 



Pn (o,T) = _n v — — / ^ ; \w ( 37 ) 

l-C-ln(A)], (38) 



where the second expression, Eq. (|3S|), is valid in the limit A < 1, A = h 2 k 2 D /k B TAm* , 
k 2 D = Aite 2 n/ ek B T is the inverse of the squared Debye screening length and C ~ 0.577 is 
the Euler's constant. Notice that, in the non-degenerate limit, p^ becomes almost inde- 
pendent of the total density n and independent of the spin density components n a , while 
a quantum mechanical dependence on H survives even in this regime, p^ tends to zero as 
(k B T)~ 3 / 2 ln(k B T) as T — > oo. The inset of Fig. |I| illustrates the comparison between p-^ and 
its asymptotic form. This approximation becomes valid for T ^> Tp, but, since Tp ~ n 2 ^ 3 , 
such limit is fulfilled only at very low carrier densities. 



D. The mixed (degenerate/nondegenerate) case 

A very interesting limit is the one corresponding to a spin polarization process, for 
which rif — > n and — > 0. This is indeed relevant for one of the problems we want to 
analyze, i.e. a semiconductor with strongly spin-polarized carriers. This system is in a very 
peculiar state: its minority down-spin population is non-degenerate, i.e. k B T ^> Epi, E Fa = 
h 2 (Qir 2 n a ) 2 / 3 /2m* , while, for low enough temperatures, the majority up-spin population is 
degenerate, i.e. k B T <C -Eft- The expression for the non-interacting spin-resolved density- 
density response functions entering the spin transresistivity Eq. (|29|) can then be taken from 
the previous sub-sections and are given by 

m* 2 uj' 

^(q,u/ ] T) = -— 3- forO<g<2fc FT (39) 
47m q 

= otherwise (40) 

and 

// / / rr* /7rm*\5 a , 3m* J 2 .u' , t , 

jfafau/lT) = - {—) exp(-^^-)-, (41) 

where Eq. (40) is valid up to first order in u' and Eq. (f4l[) represents the classical limit 
{% — > 0) of Eq. (|36|) . Using Eqs. p0|) , ( |4lD and the approximation, due to the small density 
of down-spin carriers e(q,uj') ~ ej(g, 0), the expression for the spin transresistivity becomes 

Rep ' l<,) ' T) = --—U^Tn, i FMJ/«-.o)|>- (42) 
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Eq. fl42|) is very similar to Eq. (|35|), the result obtained when both components are degen- 
erate. This is a consequence of the fact that, in the appropriate regime, the non-degenerate 
spectral function Eq. ( |4"1| ) presents the same dependence in uj' and q as the degenerate one. 
We want to underline that, in this limit, pn(0,T) is independent of nj_ and that in any case 
p^l —* for T — ► 0. 



IV. SPIN STIFFNESS OF AN INTERACTING SPIN-POLARIZED ELECTRON 

GAS 

The other ingredient entering the drift-diffusion expression for the current, Eq. is 
the longitudinal spin stiffness matrix S a p. In particular we are interested in the combination 
S = d 2 f(n,m,T)/dm 2 = (S^ — + S±± — S|f)/4, where m = n-f — which gives the 
curvature of the free energy with respect to the magnetization at constant density. This 
quantity coincides with the inverse of the longitudinal spin susceptibility of the uniform 
electron gas. 

We evaluated S numerically starting from the formulas provided by Tanaka and Ichimaru 
p6| , who calculated the free energy density of the three-dimensional electron gas as a function 



of temperature, density, and spin polarization. 

Fig.| shows S divided by its non-interacting value S ni as a function of the dimensionless 
temperature T/Tp for various densities, starting with n = 4.2 x 10 17 cm~ 3 for the upper 
curve down to n — 4.2 x 10 11 cm~ 3 for the lowest one, decreasing by an order of magnitude 
from one curve to the next. 

Two regimes are clearly visible. For densities larger than a critical value n c ~ 4.2 x 
10 13 cm -3 the spin stiffness decreases monotonically with decreasing temperature settling to 
a finite value in the ground-state (top four curves). For densities lower than n c a second- 
order ferromagnetic transition occurs: the critical temperature T c raises from ~ at n = n c 
to a sizeable fraction of the Fermi temperature at n = 4.2 x 10 11 cm -3 . As in any second- 



order transition, the spin stiffness vanishes at the transition temperature [2?|. For T < T c 
the spontaneous magnetization m is given by the stable minimum of the free energy, which 
satisfies the conditions 



df(n,m,T) 



S 



dm 
d 2 f(n,m,T) 



dm 2 



0. 

> 0. (43) 



m=m 



Obviously, ferromagnetism shows up only at extremely low densities, and the electron gas 
model may break down well before getting to such densities: it is nevertheless instructive to 
study the repercussions of the behavior of S on the spin diffusion constant both above and 
below the critical density. 

Let us now consider the limiting behaviors of the spin stiffness at high (T >> Tp) and 
low (T « Tp) temperatures. 
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A. The high temperature limit 



In the high-temperature limit (T ^> Tp(n)) the free energy density of the electron gas 
has the following expansion: 



f(n h n h T) ~ k B Tj2 n a[^(n a \^) - 1] 

a 

keT _ 2 3 

+ OS/2 2^ n a*T 



2^e 2 X 2 T Y.< 

a 

2n l/2 e 3 n 3/2 

3 {k B Ty/ 21 



(44) 



where At = {2kTi 2 /m^ksT) 1 ! 2 is the thermal wavelength. The first term is the free energy of 
the classical ideal gas, the second term is the leading quantum correction for noninteracting 
Fermions ||28|| ; the third term is the leading quantum/interaction correction, namely, the 
high temperature exchange free-energy |29j , the last term is the leading classical interaction 
correction from Debye-Huckel theory. Only the first three terms depend on the magnetiza- 
tion, and therefore contribute to the spin stiffness. Taking a second derivative with respect 
to magnetization we find, after simple calculations, 



— = 1 + — l 

Sr. n 



2 



n\ 3 T _ 8ttW/3 
2V2 k B T 



Z_ ^tn ( ^ A 3 )2/3 



(45) 



where S c is the Curie spin stiffness of an ideal classical gas of density n: 



S c = M^. (46) 

Notice that the leading interaction correction to the noninteracting spin stiffness is negative, 
in agreement with the behavior seen in Fig. 6. 



B. The low temperature limit 

Let us now examine what happens in the limit T — - > 0. Above the critical density 
n c the spin stiffness simply tends to a constant zero-temperature limit, smaller than the 
noninteracting value, in agreement with the Landau theory of Fermi liquids. 

For n < n c the low temperature phase is ferromagnetic, and in this case the density of 
majority spin electrons (n^) approaches the total density, while the density of minority spin 
electrons (n^) tends to zero for T — > 0. To understand the behavior of the spin stiffness 
shown in Fig. 6 we assume that, in the nearly 100% polarized limit the free energy can be 
written as the sum of the ground-state energy of the degenerate interacting up-spin gas plus 
the free energy of an infinitely dilute noninteracting down spin gas: 

f{n-[, ri|, T) ~ eo(^t) — ksTrii ln(l — C) + constant, (47) 
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where €o(n) is the ground-state energy density of a degenerate liquid of up-spin electrons, 
and ( = (n^ — ni)/(n-[ + ni) is the degree of spin polarization. Obviously, this approximation 
ignores the correlation between down- and up-spin electrons, or, more precisely, presumes 
that this correlation is smaller than the entropic term ksTn^ ln(l — () for £ — ► 1. 

Starting from Eq. ([IT]) it is trivial to show that the minority spin density vanishes for 
T -> as 

n { ~ ne<i {n)/kBT (48) 
where e' (n) = deo(n)/dn < at low density, while the spin stiffness goes as 



An 



1 



(49) 



which diverges exponentially for T — > 0. This is precisely what our numerical calculations, 
based on the formulas of Ref. ( ||26||), indicate. As we shall see, this result is important 
in understanding the behavior of the diffusion constant of a unipolar spin packet when the 
system is fully spin polarized. 



V. THE EVOLUTION OF A SPIN PACKET 

We now examine the motion of a spin packet under the effect of a uniform electric 
field. Let us apply Eq. (|) to calculate the time evolution of a spin packet obtained by 
injecting an excess spin density Am(f, 0) = AMi(f) near the origin at time t = 0. We 
denote by m(r, t) = n^(r, t) — 7ij(r, t) the net spin density at point r and time t, by = 
— the uniform value of the spin density at thermodynamic equilibrium, and by 
Am(f,t) = m(f, t) - m (0) the excess spin density following spin injection. To solve this 
problem we combine the expression for the current density Eq. @ with the generalized 
continuity equations for the spin-density components 

dAn a (f,t) An a (f,t) An a (f,t) -± f 

— dt — = — T t — + 77- " V ' a(r1, (50) 

where T s f >a is the spin-flip relaxation time for the a component. Substituting in this equation 
the drift-diffusion expression for the current Eq. (|]) we obtain the two equations 

dAn a (r, t) An a (f,t) An g (f, t) 

9t 7~sf,a T~sf,a 

+~ E ^V(An /3 ) + ~a a ^—- + £ [VDapVnp + D aP V 2 {An p )] , (51) 



(a =| or |) where 



P 



Eq. (|51|) includes the term a a V ■ E fe = a a (An^ +An±)/e. This term is "dangerous" because 
it contains the product of a large quantity a a times a small quantity, the space charge 
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Ari| + Ani, the product itself being of the order of the other quantities of interest in the 
calculation. 

It is tempting, but wrong, to invoke the local charge neutrality constraint 

An t (r) = -Anj(r) (53) 

at this point. Instead, we will first combine the two component of Eqs. (|5l|) to eliminate the 
space charge (V ■ E) term, and only after doing that can we impose, without serious loss of 
accuracy, the charge neutrality constraint. To eliminate the V ■ E term we multiply each 
component of Eq. (pTTD by the conductivity of the opposite channel, and, in order to get 
the equation of motion for Am, we subtract the equation for a-spins from the equation for 
a-spins. Only at this point we impose the local charge neutrality constraint Eq. ([53]). With 
this procedure we obtain the correct drift-diffusion equation for Am(f, t) 

dAm(f,t) Am(r,t) E a ^V(D Q VAm(f,i)) 
dt r s E> a 

+ Eaa f" i-VAm(f,t) (54) 

where r s = (l/r^j + l/r^) -1 is the spin relaxation time, which is very long ||4 , , 

D a = D aa — D a ai (55) 

/Iq, f^aa ^aai (56) 

E is an externally applied electric field, and the matrix [i a p is defined as 

da a/ 3 da a/ 3 

^ " ~Dn~ ~ ~Dn^ (5?) 

= ±2 a ^ , plus if a =|, minus otherwise. (58) 
dm 

This is a generalization of the familiar relation between mobility and conductivity: it takes 
into account the dependence of the mobility /j, a p on both spin density components. The 
second term in Eq. ( |57D accounts for the reduction of the mobility in the a channel due to 
the drag of the a spin population on the a population 0. 

The fact that the mobilities enter Eq. ([54]) as a spin symmetric combination (Eq. (|56f) ) 
while the diffusion constants are in a spin antisymmetric combination (Eq. (|55|)), reflects 
the fact that the electrostatic field has the same sign for both spin components, while the 
density gradients have opposite signs (see Eq. (plf)). 

If we consider the linear regime - i.e. we neglect terms of the order of (Vn a ) 2 - Eq. ([57]) 
reduces to the more familiar 

a a/3 = en a fi a/3 (59) 

and Eq. (|>|) can be written as 

dArr^tl = _ Am(f,t) + f) + j . f Am{jf f) (6Q) 

at t r 
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where 



and 



Set 0~aPa /m\ 

Us = (61) 



D s = (62) 

L-ia ®~a 



are the effective mobility and diffusion constants [|30|. We underline that in reality the range 



of validity of Eq. (p0|) extends beyond the linear approximation into the classical regime, 
i.e. to high carrier densities or high temperatures, since in that regime the relationship 
between density and conductivity is linear. Eqs ( |6~T1) and fl6~2"D show that the mobility and 
the diffusion constants of the packet are weighted averages of, respectively, the mobilities p a 
and diffusion constants D a of the two spin channels, the weight being the conductivity a a 
of the opposite channel. This is due to the local charge neutrality constraint Eq. (|53|), that 
forces the two components of the disturbance to travel together, so that the conductivity 
of each spin channel is strongly influenced by the motion of the disturbance in the other 
channel. In the non interacting limit Eqs (|6T|) and fl62| ) reduce to the expressions presented 
in §. 

The solution of Eq. (|60D is 

AmM = WD^ e ' (63) 

Eq. ( |5"5D has the form of a Gaussian packet that drifts under the effect of the electric field 
E with a pace determined by p s , and spreads in time at a rate determined by D s . The 
mobility and diffusion constants of electron-hole packets of similar shape can be measured 
through the Haynes-Shockley experiment [31]. Thus a similar experiment can in principle 



determine p s and D s independently, provided that r s is sufficiently long. 
With the help of some algebra, Eq. (Ell) can be rewritten as 



In 1 

Ps = { . 64 



where a a is given by Eq. (p2|), and Eq. ( p2|) as 

S 4 



D s = (65) 



e 2 Ea Po 

k B T S 1 



e 2 n S c [(n^ ni /n 2 )(p D1 + p Dl ) - p n ) 



(66) 



where 



Pa = Paa ~ Paa, (67) 



S c is given by Eq. fl4"6|), p Da = m*/ne 2 r a is the ordinary Drude resistivity associated with 
the a spin channel, and - a negative number - is the spin drag transresistivity discussed 



in Sec. Ill 
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Eq. (^) can be seen as the generalization to a spin packet of the ordinary relationship 
between mobility and conductivity, while Eq. (|65[) is the corresponding generalization of the 
Einstein relation. To derive Eq. (|66|), we made use explicitly of the structure of the matrix 



p a p discussed in Sec. |IIB 



The calculation of p s and D s is simplified under the assumption that the scattering times 
for the two spin components are not too different, i.e., r-j- = tj_ = Td- This assumption is well 
justified for non-degenerate carriers. In the degenerate case l/r Q has a weak n l J 3 dependence 



on the density. From Eqs. ( |6TD and (|66|) we obtain 



Vs = — 68 

m* 

and 

p s k B T S 1 

D s = ^ 1 — , 69 

e S c l- Pn/pD 

where po = m*/ne 2 TD is now the ordinary Drude resistivity. 

Eq. fl68D tells us that the mobility of the packet is not explicitly modified by Coulomb 
interaction and in fact coincides with the ordinary homogeneous mobility. Strictly speaking 
this result is only valid under the assumption that up-spin and down-spin electrons have 
equal mobilities and thus drift at the same speed in an applied electric field. Coulomb inter- 
actions, being Galilean-invariant, cannot change the total momentum of such a uniformly 
drifting electron gas. 

The situation is completely different for the diffusion constant. As the spin packet spreads 
out the up- and down-spin currents are directed in opposite directions and friction arises: for 
this reason the expression for D s contains the spin-drag resistivity as a factor that reduces the 
diffusion. In addition, the Coulomb interaction together with the Pauli exclusion principle 
reduces the energy cost of spin-density fluctuations (i.e., the spin stiffness) decreasing further 
the rate of diffusion of a spin packet. 



VI. EXPERIMENTAL OBSERVATION OF INTERACTION EFFECTS IN 

SPIN-POLARIZED TRANSPORT 

A. Direct observation of spin Coulomb drag 

We will now propose an experiment aimed at detecting the effect of the spin Coulomb 
drag and measuring the spin trans-resistivity. We describe the experiment as it could be 
done on metals, but the same scheme could be applied to semiconductors (in which the drag 
effect is larger), provided that an appropriate method of injecting spin current is used. 

The setup is shown in Fig. ||: a paramagnetic metal film of thickness L is sandwiched 
between two ferromagnets polarized in the same direction. A battery is connected to the 



ferromagnets inducing a spin-polarized current [0]] from the first ferromagnet ("injector") 
through the paramagnet and toward the second ferromagnet ("receiver"). The injector and 
receiver are chosen to be half-metals, i.e, they have only electron states of spin j at the Fermi 
level (see Fig.||). It follows that the injected current is carried only by spin j electrons 
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|32fl . If we choose L « 5 S , where S s is the spin relaxation length, we can safely neglect spin- 
flip processes. Spin relaxation lengths are relatively large in some materials (S s ~ 100/xm 
in Al ||11||), so the condition L << S s is not particularly restrictive. At any rate, we have 
seen in Section II that the coupling between spin channels induced by spin-flip scattering 
is expected to be much weaker than the one due to the Coulomb interaction. Due to the 
spin Coulomb drag, the injected will drag spin j electrons toward the junction with the 
receiver. But, since there is no conduction band available in the receiver for spin j electrons 
the circuit will behave as an open circuit for spin j electrons, i.e., j± = 0. The vanishing of 
jl is an indication that the Coulomb drag force is exactly balanced by the gradient of the 
electro-chemical potential for spin down electrons 

1 e 

= PlTJt (70) 

What Eq. ( [70]) tells us is that, due to the spin Coulomb drag, there will be a measurable 
electro-chemical potential difference eEJ = ep^j^l for spin j electrons between two points 
within the metal separated by a distance / along the direction of the current. 

To measure this potential difference a second circuit including a voltmeter of very large 
resistance is connected to the regions of the paramagnet close to the junctions (See Fig. [|). 
Our purpose is to measure E±, so this second circuit must be driven by the spin j electro- 
chemical potential only. In order to accomplish this, we propose to use as contacts two 
half-metallic ferromagnetic electrodes ("detectors"), similar to the injector and the receiver, 
but polarized in the opposite direction. In this way, for the same reasons explained before, 
the detection circuit will be "open" as far as spin f electrons are concerned, and the current 
flowing in the voltmeter will be exclusively driven by the electro-chemical potential difference 
of spin I electrons. The spin trans-resistivity will then be given by pjj_ = (AVd/I^)(A/1), 
where AVd is the voltage measured by the meter, A is the cross-section of the paramagnetic 
metal, I is the distance between the detectors, and if the current flowing between injector 
and receiver. As shown by our calculations (see Fig. ^), we expect, in metals, a resistivity 
of the order of 10 -2 fiQcm that is proportional to T 2 for ksT » uo. 



B. Haynes-Shockley experiment 



The Haynes-Shockley experiment [01] demonstrated the drift and diffusion of minority 
carriers in a doped semiconductor. The experiment allows a direct and independent measure 
of the minority carrier diffusion and mobility coefficients. After a pulse of excess carriers is 
created at some point in the semiconductor, it drifts under the action of an electric field, for 
a known distance L, after which it is monitored. By measuring the drift time and the width 
of the packet, it is then possible to compute both the mobility and the diffusion constant of 
the packet, which coincide with those of the minority carriers. 

The experiments of Ref. ( can be seen as a Haynes-Shockley-type experiment, since 
they are based on a direct monitoring in space and time of unipolar spin packets. In these 
experiments, Kikkawa and Awschalom were able to measure independently the diffusion and 
the mobility of the spin packets. In the next section we are going to compare our theoretical 
results with their experimental findings. 
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VII. PARAMAGNETIC SEMICONDUCTORS AND SPONTANEOUS 
FERROMAGNETIC TRANSITION 



We will now focus on the results of our calculations for the diffusion constant of the spin 
packet. First of all we will analyze the paramagnetic regime in which = n±. In Fig.[7] 
we plot the ratio D s /D ni , where D s represents the fully interacting calculation according 
to Eq. ( |69|) and D ni is the non-interacting diffusion constant, as in ||. The figure shows 
results for n-doped GaAs, at three different temperatures and in a range of densities that 
is relevant to the experiments of Ref. [Q. We see that the interaction correction is quite 
significant, and, in the paramagnetic regime, reduces the value of the diffusion constant as 
expected, i.e. D s /D ni < 1 always. The solid lines correspond to the calculations performed 
at a temperature T = 300K, the dashed lines to T = 20K and the dot lines to T = 1.6K. 
The curves marked with " SD" correspond to the case in which interactions in D s are taken 
into account only through the spin Coulomb drag effect. The figure shows clearly that at 
low temperatures the most important many-body contribution to the diffusion is due to the 
softening of the spin stiffness, while, already at 20K the spin drag contribution becomes 
relevant, to represent most of the interaction effects at room temperature. We see that 
D s /D n i — > 1 for high densities. This is due to the enhancement of screening in this regime, 
so that the particles tend to behave as non-interacting ones. In the high density regime 
S ~ S ni as can be seen in Fig. |5|. If the temperature is high enough, the spin drag still 
reduces the diffusion constant by a sizable amount, but eventually, even this contribution 
disappears with increasing density, and D — > D n i. 

At low density, the system enters the non-degenerate regime, so that both D —>■ D c and 
D ni — > D c , where D c = p s k B T/e is the classical non-interacting diffusion constant. This 
limit will be discussed in greater detail below. 

In Fig. [8] we plot D s /D c for the same parameters of Fig. [7|. The solid lines represent the 
fully interacting diffusion constant D s while the dashed lines represent the corresponding 
non- interacting approximation D ni . We see that, despite the significant reduction due to 
the interaction correction, D s remains still considerably larger than D c , consistent with 
experimental observations ||. 

As mentioned above, in the non-degenerate limit D/D c approaches 1. It is worthwhile to 
examine how this limit is approached. In the non-interacting theory |5| the non-degenerate 
limit is approached from above because the leading correction to D c is due to the spin 
stiffness and comes from the quantum kinetic energy, i.e. it is positive. In the interacting 
theory there is an additional exchange correction to S, which is negative, and competes with 
the quantum kinetic one (see Sec. [IV A| ) . However, the leading correction to D c does not 
come from the spin stiffness term, but from the spin Coulomb drag, and it is always negative 
(see Eq. ©). In fact, for T » T F (n), p n /p D ~ [n/{k B Tfl 2 ] \n(n/(k B T) 2 ) (see Eq. @), 
and the logarithmic term dominates over corrections entering the spin stiffness in both the 
n — > and T — > oo limits. Thus, due to interactions, D s /D c — > 1 from below always. This 
is evident in Fig.[7| where, for a fixed temperature, D s /D c becomes negative as the density 
decreases under a certain threshold, while for a fixed density the ratio becomes negative 
when shifting to curves calculated at higher temperature. 

We want to stress that D s also displays a marked dependence on the sample mobility 
that affects the diffusion constant through Drude resistivity p B . The higher the mobility, 
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the more important becomes the factor containing the spin transresistivity 1/(1 — Pu/Pd) 
(see Eq. (|S9|) ). As we already underlined in Sec. |l [I A| the ratio Ph/pd can become very 
relevant and even greater than 1 (see Fig. ^ and |3|). The diffusion constant in these cases is 
then regulated by the spin drag effect that cannot be neglected. 

Perhaps the most interesting feature of Eq. (^) is the possibility of a large variation in 
D s when the electron gas undergoes a ferromagnetic transition. From the curves in Fig.|5| 
and Fig. [5] we see that S and D s vanish at the transition temperature and increase sharply 
as the system settles in the fully polarized state. In the case of intrinsic ferromagnetism, 
the critical behavior of D s is completely due to Coulomb interactions among carriers. 

As can be seen from Fig. |9], D s /D c — ► 1 as the system fully polarizes. In this limit 
ni — > 0, Pii/pd — > 0, and S/S c — > 1, as demonstrated in Sec. [IV B| , so that D s reduces to 
the diffusion constant of carriers of minority orientation (which are non-degenerate), i.e. to 
the classical value D c = k B Tp/e. 

Unfortunately, in an ordinary electron liquid, the ferromagnetic transition is predicted 
to occur only at extremely low densities. There is, however, an interesting variant: semi- 
conductors doped with magnetic impurities (for example Mn) can undergo a ferromagnetic 
transition at rather high carrier densities 0,0, n ~ 10 20 cm -3 for (Ga,Mn)As, and temper- 
atures T < T c ~ 110K 33] . Our theory on the dynamics of a spin packet can be extended 



to these systems with similar results |3q|. This extension will not be pursued here. 



VIII. CONCLUSIONS 



In this paper, we have tried to demonstrate the importance of many-body effects in spin 
polarized transport. We have discussed in detail the spin Coulomb drag effect, an intrinsic 
source of friction in spin transport that can limit spin currents even in the purest materials. 
We have worked out the behavior of the spin transresistivity pjj in different physical regimes 
and shown that it ranges from 10~ 8 Qcm in metals to 10~ 3 — 10~ 2 Qcm in semiconductors. 
Moreover, the ratio Ph/pd, which is only a few percents in metals, becomes comparable 
to, or even larger than unity in semiconductors. We hope that an experimental group will 
soon take up the challenge of measuring the spin transresistivity, for example through the 
experiment we suggest, in order to confirm the theory. 

We have also demonstrated the importance of including Coulomb interactions in a quan- 
titative theory of spin diffusion, and shown that a measure of D s for a unipolar spin packet 
would be a sensitive probe of many-body effects such as the spin-Coulomb drag and the 
Coulomb enhancement of the spin susceptibility. 

Finally we have studied the behavior of D s at and below a spontaneous ferromagnetic 
ordering transition and found that D s exhibits a critical behavior. 
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FIGURES 



FIG. 1. Spin transresistivity p-M as a function of temperature (rescaled by Tp) for GaAs param- 
eters (m* = 0.067m e , e = 12, /u = 3 x 10 3 cm 2 /Vs). Each curve corresponds to a different density: 
n\ = 1.5 x 10 16 cm -3 , 7i2 = 1.5 x 10 17 cm -3 = 1.5 x 10 18 cm -3 . Inset: comparison between p+\ 
and its analytical approximation in the non-degenerate regime vs temperature (rescaled by Tp) for 

n = 722. 

FIG. 2. Upper panel: Ph/pd as a function of temperature for GaAs parameters (m* = 0.067m e , 
e = 12, n = 1.5 x 10 17 cm~ 3 ). Each curve corresponds to a different mobility: A = W 2 cm 2 /Vs, 
B = 10 3 cm 2 /Vs, C = 3 x 10 3 cm 2 /Vs D = 10 4 cm 2 /Vs, as labelled. Lower panel: same as upper 
panel but for (Ga,Mn)As parameters (m* = 0.5m e , e = 12, n = 1.2 x 10 19 ). 

FIG. 3. P\\l Pd as a function of carrier density for GaAs parameters (m* = 0.067m e , 
e = 12). The solid curves are calculated at the temperature T = 300-ftT, while the dashed curves at 
T = 20K. For each temperature two different mobilities are considered, C = 3 x W 3 cm 2 /Vs and 
D = 10 4 cm 2 /Us, as labelled. 

FIG. 4. Temperature and density dependence of p-t\ in the degenerate regime. The top line 
corresponds to r s = 7. The electron-gas parameter is decremented by 1 starting from the top. 

FIG. 5. Spin stiffness S, rescaled by its non-interacting approximation S n i vs T /Tp. The carrier 
density is n = 4.2 x 10 11 cm -3 for the lower curve and increases by a factor 10 for each line starting 
from the bottom. The cusps correspond to the onset of ferromagnetism. 

FIG. 6. (a) Experimental setup to detect the spin Coulomb drag effect: the voltage AU is 
applied between two parallel half-metallic ferromagnets (injector (inj.) and receiver (rec.)) that 
sandwich a paramagnet (P). The voltage AVd is detected using two ferromagnetic electrodes (d) 
similar to the injector and the receiver, but polarized in the opposite direction, (b) Schematic 
band-structure of injector, receiver, detectors and paramagnet P. (c) Schematic behavior of the 
chemical (£|) and electro-chemical potentials ip a (f) (see Eq. ||). For this setup £j is basically 
constant (not shown in the figure). 

FIG. 7. The interacting diffusion constant of a spin-packet D s rescaled by its non-interacting 
approximation D n i vs density for different temperatures: solid lines correspond to T = 300i^, 
dashed lines to T = 2QK and dotted lines to T = 1.6K. For each temperature, we plot also the 
curve obtained considering interactions only through the spin Coulomb drag effect (labeled as SD). 
In all the calculations the dielectric constant of the semiconductor is e = 12 and the mobility is 
H = 3 x 10 3 cm 2 /Us. 
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FIG. 8. The interacting (solid line) and non- interacting (dashed line) diffusion constant of 
a spin-packet rescaled by its classical non-interacting approximation D c vs density for different 
temperatures. In all the calculations the dielectric constant of the semiconductor is e = 12 and the 
mobility is \i = 3 x 10 3 cm 2 /Vs. 

FIG. 9. The diffusion constant of a spin packet rescaled by its classical non-interacting ap- 
proximation D c vs temperature in a low-density electron gas with no magnetic impurities. 
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